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Phase structure of the (2+l)-dimensional model with four-fermion interaction of spin-1/2 quasi- 
particles (electrons) both in the fermion-antifermion (or chiral) and fermion-fermion (or super- 
conducting) channels is considered at nonzero chemical potential [s, and under the influence of an 
in-plane, i.e. parallel to a system sheet, external magnetic field B\\ . It is shown that at sufficiently 
large values of p and/or B« the Cooper pairing (or superconducting) phase appears in the system at 
_ arbitrary relation between coupling constants, provided that there is an (arbitrary small) attractive 

interaction in the superconducting channel. In particular, at sufficiently weak attractive interaction 
in the chiral channel, the Cooper pairing occurs even at infinitesimal values of fj, and/or Bp The 
£N) _ superconducting phase of the model is always a paramagnetic one. 
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INTRODUCTION 



Last years much attention has been paid to investigation of (2+l)-dimensional quantum field theories (QFT) and, 
in particular, to models with four-fermion interactions of the Gross-Neveu (GN) [l| type. Partially, this interest is 
explained by more simple structure of QFT in two-, rather than in three spatial dimensions. As a result, it is much 
easier to investigate qualitatively such real physical phenomena as dynamical symmetry breaking Ifi-lp and color 
£Lj superconductivity [lOj as well as to model phase diagrams of real quantum chromodynamics (QCD) etc. in the 
framework of (2+l)-dimensional QFT. Another example of this kind is the spontaneous chiral symmetry breaking 
induced by external magnetic fields. This effect was for the first time studied also in terms of (2+l)-dimensional 
GN model [T2J, |l3j. Moreover, these theories are very useful in developing new QFT techniques like the optimized 
perturbation theory [lj [l4j , and so on. 

However, there is yet another more serious motivation for studying (2+l)-dimensional QFT. It is supported by the 
fact that there are many condensed matter systems which, firstly, have a (quasi-)planar structure and, secondly, their 
excitation spectrum is described adequately by relativistic Dirac-like equation rather than by Schrodinger one. Among 
these systems are the high-T c cuprate and iron superconductors [15j |. the one-atom thick layer of carbon atoms, or 
graphene, [l6l . [l7j etc. Thus, many properties of such condensed matter systems can be explained in the framework 
of various (2+l)-dimensional QFT, including the GN-type models (see, e.g., p^| - |26j and references therein). 

In the recent papers [24T-|26| a competition between chiral symmetry breaking (excitonic pairing) and supercon- 
ductivity phenomenon (Cooper pairing) was investigated in the framework of (2+l)-dimensional GN-type models. 
There the influence of such external factors, as temperature T, chemical potential /j, and external magnetic field B± 
perpendicular to the system plane, on the chiral and electromagnetic U(l) symmetries was studied. In particular, it 
was shown in 25] that sufficiently strong perpendicular magnetic field B± destroys the superconducting state of a 
planar system. 

In the present paper the (2+l)-dimensional GN-type model, which describes four-fermion interaction of quasipar- 
ticles with spin 1/2 (electrons) both in the chiral (with coupling constant G\) and Cooper pairing (with coupling 
constant G2) channels, is considered at T — and // / 0. In addition, we suppose that the planar system of electrons 
is subjected to an external magnetic field B = B± + B\\ in such a way that B± = (B\\ is the in-plane, i.e. parallel to 
the system plane, magnetic field). As a result, in our consideration the external magnetic field B couples only to the 
spin of electrons and not to its orbital angular momentum. (The interaction between B and orbital angular momentum 
of electrons appears in planar systems only at B± 7^ 0.) In this case the way to account for the influence of external 
parallel magnetic field on the phase structure of any planar system is to introduce the Zeeman terms into Lagrangian, 
which are really due to Zeeman effect. These terms modify effectively the chemical potentials corresponding to the 
electrons with different spin projections along the direction of an external magnetic field. I n doing so we show that 
the external magnetic field -By parallel to the system plane induces the superconductivity phenomenon in the initial 
model even at infinitesimal coupling Gi of the Cooper pairing channel. If in addition the coupling constant G\ is 
sufficiently small, then superconductivity appears at arbitrary weak in-plane magnetic field B\\ (of course, provided 
that G 2 > 0). 



In the recent paper 12 II the magnetization of planar systems exposed to an external in-plane magnetic field was investigated in the 
framework of the same (2+l)-dimensional GN model but in the particular case Gi 7^ 0, G2 = 0. In particular, it was shown there that 
at sufficiently high values of By there is a restoration of the chiral symmetry of the model. 
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The paper is organized as follows. In Sec. II the (2+l)-dimensional GN-type model with four-fermion interactions 
in the fermion-antifermion (or chiral) and fermion-fermion (or superconducting) channels is presented. Here the 
unrenormalized thermodynamic potential (TDP) of the model is obtained in the mean field approximation, taking 
into account the nonzero values of the chemical potential \x and the external in-plane magnetic field B» (temperature 
is put equal to zero). In the next Sec. Ill a renormalization group invariant expression for the TDP is obtained 
whose global minimum point provides us with chiral and Cooper pairing condensates. In Sec. IV phase structure of 
the model is described at /i = and B» ^ 0. In particular, it is found out in this section that infinitesimal values 
of -By induce the superconductivity phenomenon in the case of a rather weak attractive interaction in the fermion- 
antifermion channel. Finally, in Sec. V the (/^, |-By |)-phase diagrams are presented for some representative values of 
coupling constants. We show in this section that for arbitrary relations between coupling constants superconductivity 
is induced in the system at sufficiently large values of /i and/or in-plane magnetic field B». Some related problems of 
our consideration are relegated to two Appendices. 



II. THE MODEL AND ITS THERMODYNAMIC POTENTIAL 

We investigate the influence of an external magnetic field B on the (2+l)-dimensional GN-type model with 
two massless fermions ipx (x) and ip2(x) belonging to a reducible 4-component spinor representation of the (2+1)- 
dimensional Lorentz group (the spinor fields tpi{x) and ip2{%) are introduced for quasiparticles (electrons) with spin 
projections 1/2 and -1/2 on the direction of the magnetic field B, respectively). The model describes low-energy dy- 
namics of quasiparticles both in the fermion-antifermion (or chiral) and fermion-fermion (or Cooper pairing) channels. 
If external magnetic field is parallel to the system plane, i.e. B = B\\ , then the Lagrangian has the following form 

2 

L = E 

k=X 





where \x is a fermion number chemical potential and the ^-term is introduced in order to take into account the Zeeman 
interaction energy of electrons with external magnetic field -By. Hence, in our case v = g[iBB/2, where B = \B»\, 
g is the spectroscopic Lande factor (below it is supposed throughout the paper that g = 2) and \ib is the Bohr 
magneton. In (1) the symbol T denotes the transposition operation and the summation over the repeated Lorentz 
indices p = 0,1,2 is implied. Moreover, C = 7 2 is the charge conjugation matrix. The algebra of the 7 p -matrices 
as well as their particular representations are given in Appendix [X] The model (1) is invariant under the discrete 
chiral transformation, ipk — > 7 5 V'fc (the- particular realization of the 7 5 -matrix is presented in Appendix |A|) . as well 
as with respect to the transformations from the continuous U(l) fermion number group, ipf. —> exp(ia)i/jk {k = 1,2), 
responsible for the fermion number conservation or, equivalently, for the electric charge conservation law in the system 
under consideration. 

The linearized version of Lagrangian ([1]) that contains auxiliary bosonic fields <j{x), A(:r), A* (a;) has the form 



a 1 A* A v - 

C = h > 

AGx 4G 2 f-f 



2 



(2) 



where fix = A* + v, /j,2 = fJ, — v and from now on v = (IbB. Clearly, the Lagrangians ([T]) and @ are equivalent, as can 
be seen by using the Euler-Lagrange equations of motion for scalar bosonic fields which take the form 
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<j{x) = -2G!^Vi fe Vfc, A(z) = -2<7 2 ^V>*C%, A*(x) = -2G 2 J2^kC^. (3) 

fc=l k=l k=X 

One can easily see from Q that the neutral field cr(x) is a real quantity, i.e. (cr(.x))^ = a(x) (the superscript symbol 
f denotes the Hermitian conjugation), but the (charged) difermion fields A(x) and A* (a;) are mutually Hermitian 
conjugated complex quantities, so (A(x))^ = A*(x) and vice versa. Moreover, with respect to parity transformation 
P (see also the comment in Appendix A), 

P : ip k (t,x,y) ^ 7 5 7Vfc(t,-z,y), fc = l,2, (4) 

the fields tr(a;), A(x) and A*(x) are even quantities, i.e. they are scalars. If the difermion field A(a;) has a nonzero 
ground state expectation value, i.e. (A(x)) =^ 0, the Abelian fermion number £7(1) symmetry of the model is 
spontaneously broken down and the superconducting phase is realized in the model. However, if (cr(x)) ^ then the 
discrete chiral symmetry of the model is spontaneously broken. 

Let us now study the phase structure of the four-fermion model (1) by starting from the equivalent semi-bosonized 
Lagrangian In the mean-field or, equivalently, in the one-loop approximation, the effective action S c ff (cr, A, A*) 
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of the considered model is expressed by means of the path integral over fermion holds: 
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exp(z5 cff (a,A,A*)) = f JJ[#j][#,] exp(i f £d 3 x), 
^ 1=1 J 



where 



5 off (a,A,A*) 



r a 2 {x) A(x)A*(x) 



d A x 



cir- 



4Gi 4G 2 

The fermion contribution to the effective action, i.e. the term S B s in ([5]), is given by: 

A 



(5) 



exp(iS off ) = / exp|z / ^ ipk{l P id l 

J 1=1 k=l L 



J x}. (6) 



The ground state expectation values (<r(x)), (A(x)), and (A*(x)) of the composite bosonic fields are determined by 
the saddle point equations, 



5S cS 



= 0, 



£5, 



off 



= 0, 



= 0. 



(7) 



5a(x) ' 5A(x) ' SA*(x 

For simplicity, throughout the paper we suppose that the above mentioned ground state expectation values do not 
depend on spacetime coordinates, i.e. 

(a(x))=M, (A(x))=A, (A*(x))=A*, (8) 

where M, A, A* are constant quantities. In fact, they are coordinates of the global minimum point of the thermo- 
dynamic potential (TDP) fi(M, A, A*). In the mean-field approximation this quantity is defined by the following 
expression: 



J d 3 xtt{M, A, A*) = -S eS {a(x), A(x), A* (x)} 



<7(x)=A?r,A(x)=A,A*(x)=A* 



which gives 



d 3 xft(M, A, A* 



d 6 x 



M 2 AA* 



4Gi 4G 2 



+ i 



In I / JJ[#;][#;]exp U I ^ \^k D hi>k 

\^ (=1 J k=l 



d A x 



(9) 



where Dk = 7 iy «<9„ + (J^kJ — M ■ To proceed, let us first point out that without loss of generality the quantities A, A* 
might be considered as real ones. So, in the following we will suppose that A = A* = A, where A is already a real 
quantity. Then, in order to find a convenient expression for the TDP it is necessary to invoke Appendix B, where 
the path integral similar to (|9]) is evaluated. So, taking into account in ([9]) the relation (|B7I) we obtain the following 
expression for the TDP of the GN model (1) at zero temperature, T = 0: 



M 2 A 2 
W A ^4GT + 4G^ 



E 



d 3 p 

(2^ 



In 



(10) 



where (£ | k ) 2 = E 2 + fi 2 + A 2 ± 2^M 2 A 2 + fi 2 E 2 and E = M 2 + \p\ 2 . Throughout the paper we suppose that 
A* > 0, v > 0, M > and A > 0. Using in the expression (fTUf a rather general formula 



dp ln (po -A) = m\A\, 



where A is a real quantity, it is possible to reduce it to the following one: 



fl(M, A) ee Q un (M, A) 



M 2 A 2 
4GT + AG 2 



2 * ^2 

Ey j2^{ £ t k + £ A,, 

k—l 



(12) 



The integral term in (|12[) is ultraviolet divergent, hence to obtain any information from this expression we need to 
renormalize it. 



2 Otherwise, phases of the complex values A, A* might be eliminated by an appropriate transformation of fermion fields in the path 
integral JS}. 
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III. THE RENORMALIZATION PROCEDURE 



First of all, let us regularize the zero temperature TDP (| 12[) by cutting momenta, i.e. we suppose that \p±\ < A, 
\p2\ < A in (H2J). As a result we have the following regularized expression (which is finite at finite values of A): 



n re9 (M,A) 



M 2 A 2 _ _1_ A 



k=l 



dpi 



dp 2 {£ 



'A, ft 



£ 



A, ft 



Let us use in (IT31 the following asymptotic expansion (k = 1,2) 

M 2 + A 2 



°A,k 



£ 



A, ft 



0(l/\p\% 



(13) 



(14) 



where \p\ = \Jp\ + p\- (Note, the leading asymptotic terms in (| 14[) do not depend on /ii,2-) Then, upon integration 
there term-by-term, it is possible to find 



n reg (M,A) = M 2 



1 4Aln(l + V2) 
AG~i ~ 2 

1 4Aln(l + v / 2) 
AG~ 2 V 2 



4A 3 (V2 + ln(l + s/2)) 



3tt 2 



■O(A ) 



(15) 



where O(A ) denotes an expression which is finite in the limit A — > oo. Second, we suppose that the bare coupling 
constants G\ and Gi depend on the cutoff parameter A in such a way that in the limit A —¥ oo one obtains a finite 
expression in the square brackets of (|15p . Clearly, to fulfil this requirement it is sufficient to require that 



1 

4G7 



1 



4Gi(A) 



4Aln(l + V2) 1 
it 2 ngi 



1 



1 



4Aln(l + V2) 



4G 2 (A) 



1 



(16) 



where are finite and A-independent model parameters with dimensionality of inverse mass. Moreover, since bare 
couplings G\ and G2 do not depend on a normalization point, the same property is also valid for 31.2- Hence, taking 
into account in (|13l) and (fTSj) the relations (ITBl and ignoring there an infinite M- and A-independent constant, one 
obtains the following renormalized, i.e. finite, expression for the TDP 



n ren (M, A) 



lim 

A— >oo 



n reg (M 7 A) 



4A 3 (72 + ln(l + y2)) 



Gi=Gi(A),G 2 =G 2 (A) 



37T 2 



(17) 



It should also be mentioned that the TDP (|17j) is a renormalization group invariant quantity. 

Suppose that /x = and 2/ = /is-B = 0. In this case ^1,2 = 0, so the O(A ) term in (fT5"]) can be calculated explicitly. 
As a result, we have for the TDP in this particular case the following expression 



V(M, A) = n ren {M,A) 



^=0,^=0 



M 2 



A 2 (M + A) 3 

TTg 2 37T 



\M 



3?r 



(18) 



Now, let us obtain an alternative expression for the renormalized TDP (1171) at /1 7^ and v 7^ 0, i.e. at /11.2 
/x± v 5^0. For this purpose one can rewrite the unrenormalized TDP Q, m (M, A) (JT2J) in the following way 



n un (M,A) 



M 2 

2 

E 

fc=l 



A 2 
4G 2 



d 2 p 



^ 7 (2tt) 2 I A . fc W=o 



k=l 



y A.k 



£ 



A, ft 1^=0 



' £ a 1. £ 



+ 



-£: 



A, ft A,fc| Mi „ = o A, fe 1^,^=0 



(19) 



where for each k = 1, 2 



A, ft I f_i,u=o + ^A,ft I ftI y=o 



VN 2 + (Af + A) 2 + VIpI 2 + (M - A) 2 



Since the leading terms of the asymptotic expansion (|14[) do not depend on fii,2, it is clear that the integrals in the 
last sum in (|19j) are convergent. Other terms in ()19|) form the unrenormalized TDP of the particular case with /1 = 
and v — which is reduced after renormalization procedure to the expression (|18p . Hence, after renormalization we 
obtain from (|T9]) the following finite expression (evidently, it coincides with renormalized TDP (|17p): 



fT en (M,A) = V(M 



' A )~/ 70 E + £ A,k VW + (M + A) 2 - VbT + (M - A) 2 ) , 



(20) 
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III - SUPERCONDUCTIVITY 
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FIG. 1. The (gi, <72)-phase portrait of the model at 
/i = and B — 0. The shorthands I, II and III de- 
note the symmetric, the chiral symmetry breaking and 
the superconducting phases, respectively. In the phase II 
(a) — —1/gi. In the phase III (A) = —\jg%. On the curve 
L= {(gi,g2) : gi = 32}, where #1,2 < 0, the TDP minima 
corresponding to the phase II and III are equivalent. 



/Jslsil-B 




52/|ffl| 



FIG. 2. The case n = 0: The ((72, B)-phase portrait of 
the model and critical in-plane magnetic field B c (g2) at 
arbitrary fixed g\ < 0. On the critical curve B e (S2) a first 
order phase transition occurs from the chiral symmetry 
breaking phase II to the superconducting phase III. 



where V(M, A) is presented in (fT8|) . The integral terms in (|20p can be explicitly calculated. As a result, we have 



n«« ( M,A) = ^ + — + e{ t- ( m 

Kgi ng 2 \ 



- — t 



M - 



4tt^ 

(m| ~ M 2 )A 2 
47r|/i fe | 



1 

47T fe 



M • 



In 



t+ + \fi k \(M +V /M 2 fe + A 2 ) 



(21) 



where = M + A 2 ± It is not so evident, but at /i fe = (k = 1,2) the expression (JUl) for fi re "(A/, A) 
coincides with V(M, A) (JTSJ. 



IV. PHASE STRUCTURE OF THE MODEL AT fi = 



As was mentioned above, the coordinates of the global minimum point (Mo, Ao) of the TDP f2 rerl (M, A) define 
the ground state expectation values of auxiliary fields <j(x) and A(x). Namely, M = {cr{x)) and A = (A(x)). 
The quantities Mo and Ao are usually called order parameters, or gaps, because they are responsible for the phase 
structure of the model or, in other words, for the properties of the model ground state (see also the comment after 
(HI). Moreover, the gap Mo is equal to the dynamical mass of one-fermionic excitations of the ground state. As a 
rule, gaps depend on model parameters as well as on various external factors. In our consideration the gaps Mo and 
Ao are certain functions of the free model parameters 171 and <?2 and such external factors as chemical potential fj, and 
external in-plane magnetic field B. 



A. The case B = 



First of all, let us discuss the phase structure of the model (1) in the simplest case when fi — and B = 0. The 
corresponding TDP is given in (I18[) by the function V(M, A). Since the global minimum of this function was already 
investigated in [27:], although in the framework of another (2+l)-dimensional GN model, we present at once the phase 
structure of the initial model (1) at [i = (see Fig. 1). 

In Fig. 1 the phase portrait of the model is depicted depending on the values of the free model parameters g\ and 
52- There the plane (51,32) is divided into several areas. In each area one of the phases I, II or III is implemented. 
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FIG. 3. The case /i = 0: Superconducting gap Ao = 
A C rit(S2) vs (72 which is generated at critical magnetic 
field, i.e. at B — B c (g2), at arbitrary fixed g\ < 0. 
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FIG. 4. TTie case fi = 0: The magnetization m = m c (<j2) 
vs (72 which is generated at critical magnetic field B c (g2) 
at arbitrary fixed jji < 0. At B < B c (g2) magnetization 
m is equal to zero. 
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FIG. 5. 27ie case fi — 0: Superconducting gap Ao and 
magnetization m vs B at arbitrary fixed gi < and 52 = 
0.5 1 (71 1 . Curves 1 and 2 are the plots of the dimensionless 
quantities jgi|Ao and \gi\ 2 m/2fiB, correspondingly. Here 
fj,B\gi\B c = fi B \gi\B c (g 2 = 0.5| S i|) w 0.995 and | fll |A c = 
|ffi|A c ( fl2 =0.5| ff i|) « 0.098. 



FIG. 6. The case fi = 0: Superconducting gap Ao and 
magnetization m vs B at arbitrary fixed gi < and 
g-i = — 1.5|<7i|. Curves 1 and 2 are the plots of the dimen- 
sionless quantities |gi|Ao and \gi\ 2 rn/2nB, respectively. 
Here ^ s |gi|_B c = fi B \gi \B c (g 2 = — 1.5|pi |) w 0.545 and 
| 5 i|A c = |.gi|A c ( 52 = -1.5| 5 i|) w 0.838. 



In the phase I, i.e. at gi > and 32 > 0, the global minimum of the effective potential V(M, A) is arranged at the 
origin. So in this case we have Mq — (cr(x)) = and Ao = (A(x)) — 0. As a result, in the phase I both discrete chiral 
and continuous electromagnetic U(l) symmetries remain intact and fermions are massless. Due to this reason the 
phase I is called symmetric. In the phase II, which is allowed only for g± < 0, at the global minimum point (Mo, Ao) 
the relations Mq = —1/gi and Ao = are valid. So in this phase chiral symmetry is spontaneously broken down 
and fermions acquire dynamically the mass Mo. Finally, in the superconducting phase III, where 172 < 0, we have the 
following values for the gaps M = and A = — 1/<?2- 

Note also that if g\ = g% = g and, in addition, g < (it is just the line L in Fig. 1), then the effective potential 
(fT8| has two equivalent global minima. The first one, the point (Mo = — 1/g, Ao = 0), corresponds to a phase with 
chiral symmetry breaking. The second one, i.e. the point (Mo = 0, Ao = —1/g), corresponds to superconductivity. 

Clearly, if the cutoff parameter A is fixed, then the phase structure of the model can be described in terms of bare 
coupling constants Gi,G2 instead of finite quantities g\,gi- Indeed, let us first introduce a critical value of the bare 
couplings, G c — lgA ln ^ 1+v /2) • Then, as it follows from Fig. 1 and ([To]) , at G\ < G c and Gi < G c the symmetric 

phase I of the model is located. If G\ > G c , G2 < G c (G± < G c , G2 > G c ), then the chiral symmetry broken phase 
II (the superconducting phase III) is realized. Finally, let us suppose that both G\ > G c and G2 > G c . In this case 
at Gi > G2 (Gi < G2) we have again the chiral symmetry broken phase II (the superconducting phase III). 
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B. The case B^O 



In this section we study the influence of nonzero external magnetic field B parallel to the system plane on the 
phase structure of the model (1) at fi = Q. Numerical and analytical investigations of the TDP (|2"Tj) in this case show 
that its minimum points are of the form (M ^ 0, A = 0), (M = 0,A^0) or (M = 0, A = 0) only. So to study the 
properties of the global minimum point (M , A ) of the function (f2Tj) it is enough to consider its reductions on the 
M- and A-axes, where the TDP (f2"Tj) becomes (recall that v = HbB) 



6nQ ren (M, A) 



6irQ ren (M, A) 



A=0,p=0 



= 6ttcji(M) = 



M=0,/i=0 



67ru;2(A) 




(22) 



(23) 

respectively. Comparing the minima of the functions (|22[) and (|23j) . it is possible to find the global minimum point 
of the whole TDP (|21[) and its dependence on the model parameters B, g\ and g 2 , i.e. to determine the phase 
structure of the model in the particular case fi = 0. In addition, in the present Section we will study the behavior of 
a magnetization m at /i = in different phases when B varies, 

dfl ren (M,A) 



dB 



(24) 



M=M ,A=A o ,A'=0 

Since the global minimum point (GMP) of the TDP (J2TJ) coincides with the GMP either of the function wi(M) J22]) 
or the function oj 2 (A) (j2"3")h it is clear that in the chirally broken phase Ao = and the gap Mo does not depend on 
the parameter <?2- Correspondingly, in the superconducting phase we have Mq = and the gap Ao does not depend 
on the parameter g\. So, one can use the following expressions for the magnetization m in the chiral symmetry broken 
II and superconducting III phases: 

dwi(M) 

M=M 



I phase II 



= —/-IB 



m 



phase III 



= -Ms - 



dv 



= M(^_Af*)0(i/-M o ), 



(25) 



V-B 



A = An 



v\ v z + 



Agin- 



+ y^ 2 + A^ 
A 



(26) 



where 8(x) is the Heaviside step-function. 

The case g\ < 0. First of all, let us suppose that g\ is fixed and negative, i.e. g\ < 0. In this case there is a strong 
attractive interaction in the chiral channel, G\ > G c . Then it is easy to show that for arbitrary value of g 2 there exists 
a critical (in-plane) external magnetic field B c {g 2 ) (see Fig. 2) such that at B < B c {g2) the system is in the chiral 
symmetry breaking phase II (if B c (g 2 ) > 0), and it is in the superconducting phase III at B > B c (g 2 ). In other words, 
if B < B c (g 2 ) 7^ 0, then the global minimum of the TDP lies at the point (Mo = — Ao = 0) which does not 
depend on B in the interval < B < B c (g 2 ). However, at B = B c (g2) it jumps to the point (Mo = 0, Ao = A c (g2)), 
where A c (<72) vs 172 is depicted in Fig. 3. Hence, at the critical value B = B c (g2) a first order phase transition occurs 
and a superconducting gap Ao = /S. c (g2) is dynamically generated. It turns out that Ao vs B is an increasing function 
in the interval B > B c (g-i). In particular, the behavior Ao vs B is presented in Fig. 5 (at g 2 = 0.5|<7i|), Fig. 6 (at 
g 2 = — 1.5|gi|) and Fig. 7 (at g 2 = — 0.5|gi|) as the curve 1. 

Moreover, it is clear from Fig. 2 that at B < B c (g 2 ), i.e. in the phase II, the magnetization m is equal to zero. To 
explain this circumstance, recall that in the phase II the gap Mo is equal to l/|gi|. So, for all (72-values the relation 
HBB c (g 2 ) < M is valid (see Fig. 2). As a result, throughout the phase II, where v = \xbB < HBB c (g 2 ), we have 
v < Mq and hence, as it follows from the relation (|251) . the zero magnetization, mu — 0. However, when B reaches 
its critical value, B — B c {g 2 ), the superconducting phase III is generated and, as a result, the nonzero magnetization 
^0(52) appears dynamically in the system (see Fig. 4). Further growth of external magnetic field is accompanied 
by increase of the magnetization m vs B. (Evidently, in this case magnetization m must be calculated with the 
help of the expression (f2"6"|) .) For example, in Figs 5-7 at the same representative relations between gi and g 2 the 
magnetization m vs B is depicted as a monotonically increasing curve 2. 

The case g\ > 0. Now the model phase structure consideration for a positive gi-values is in order. Recall, in this 
case we have a rather weak attractive interaction in the chiral channel, i.e. G\ < G c . Evidently, if in addition 172 < 0, 



3 All the Figs. 2-10 are drawn in terms of dimensionless quantities which are obtained after multiplication of appropriate powers of |cji| 
with corresponding dimensional quantities. For example, there instead of fi, fisB, Ao, 52 we use their dimensionless analogies 
^b\9i\B, |gi|Ao, S2/|9l|- Instead of magnetization m the dimensionless quantity gfm/^ig is depicted there etc. 
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0.5 1 1.5 2 



FIG. 7. The case /i = 0: Superconducting gap Ao and 
magnetization m vs B at arbitrary fixed gi (both at gi < 
and gi > 0) as well as at 172 = — 0.5|gi|. Curves 1 and 
2 are the plots of the dimensionless quantities |<?i|Aq and 
3\gi\ 2 m/ 'hb, respectively. 



(i.i- 




0.5 1 1.5 2 



FIG. 8. The case fi = 0: Superconducting gap Ao and 
magnetization m vs B at arbitrary fixed g\ > as well as 
at <72 = 0.5<?i. Curves 1 and 2 are the plots of the dimen- 
sionless quantities |gi|Ao and \gi\ 2 m/2/iB, respectively. 



then in this case the superconducting phase is realized even at B = (see Fig. 1) due to a rather strong interaction 
in the fermion-fermion channel (G2 > G c ). Our investigations of the TDP show that in this case the superconducting 
gap Ao is enhanced by nonzero values of external magnetic field B > 0. The situation is illustrated by Fig. 7, where 
Ao and magnetization m as increasing functions of B are depicted in the particular case gi = — 0.5<7i for g\ > 0. 
Moreover, as it is clear from Fig. 7, the same behavior for Ao and m vs B remains valid for the case gi = — . 5 1 <7i | 
and negative values of g\. To explain this fact, it is necessary to take into account the remark made after formula 
(|24j) that the superconducting gap does not depend on the coupling g\ but only on 172 ■ So, it is no wonder that the 
plots of Ao and m are not changed when the parameter g\ changes the sign. 

Recall, if g% > and gi > 0, then we have at B = the phase I without any symmetry breaking, where the gaps 
Ao and Mo vanish (see Fig. 1). However, our analysis shows that at arbitrary small nonzero B the global minimum 
point of the TDP (l2lT) moves from (Mo = 0, Ao = 0) to the point (Mo = 0, Ao 7^ 0). Hence, at positive values of g\ 
and gi a continuous second order phase transition occurs from symmetric phase I to superconducting phase III when 
in-plane external magnetic field B acquires an arbitrary small nonzero value. The typical behavior of the gap Ao and 
magnetization m vs B in this superconductivity region is depicted in Fig. 8. 

In summary, we can say that at T = and [i = in the framework of the model (1) external in-plane magnetic 
field B promotes superconductivity for arbitrary relations between coupling constants g±_i (or, equivalently, Gi,2)- As 
this takes place, when at B = superconductivity is caused due to a rather strong interaction in the fermion-fermion 
channel (Gi > G c or gi < 0), the external in-plane magnetic field B > enhances this phenomenon, i.e. makes the 
gap Ao larger. However, if at B = superconductivity is absent, there is a critical value B c of in-plane magnetic 
field at which superconductivity is generated. In particular, at g± 2 > (or G1.2 < G c ) the critical field B c is an 
infinitesimal quantity. 



V. PHASE STRUCTURE OF THE MODEL ATfi^O AND B / 

The starting point of our investigations in this case is the TDP (|21j) . The behavior of the global minimum point 
(Mo, Ao) of this TDP vs fi and B supplies us with the phase structure of the model. Moreover, we are interested in 
considering such thermodynamic quantities as particle density n, magnetization m and magnetic susceptibility x> 

cW e "(M ,A ) dfr e "(M ,A ) dm 
n = ' W = dB ' X =dB- (2?) 

In the framework of the model (1) these quantities can be presented in the following form 

n = ni+ni, m — — n 2 ), (28) 

where 

dSl rm (M ,A ) drr e "(Mo,A ) . ^ e "(M ,A ) 

n\ = , n 2 = = -signal) — : (29) 

<9Mi 9^2 d\ni\ 
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FIG. 9. The case fi 7^ 0; Superconducting gap Ao and 
particle density n vs B at arbitrary fixed gi > as well 
as at g2 — 0.5<?i and /1 = 0.5/gi. Curves 1 and 2 are 
the plots of the dimensionless quantities giAo and g\n, 
respectively. 



FIG. 10. The case fi 7^ 0: Magnetization m and mag- 
netic susceptibility \ vs B at arbitrary fixed gi > as 
well as at 172 = 0.5gi and /1 = 0.5/gi. Curves 1 and 2 
are the plots of the dimensionless quantities gfrn/fXB and 
giX/A, respectively. 



are densities of particles with spin projection 1/2 and —1/2, respectively, and sign{x) is the sign-function. Note, in 
the particular case with ^ / and B = we have ni = n>2- Therefore, in this case n 7^ 0, 771 = 0. However, in the 
opposite particular case with /i = and B 7^ the relations n.1 = —772 and, as a result, 77 = 0, m 7^ are valid. It 
follows from these formulae that 



X 



2 

Mb 



d 2 O ren (M , A ) <9 2 fr e "(M ,A ) 



(30) 



(3 M i) 2 (5| M2 |) 2 

As in the particular case with /i = 0, B 7^ (see the previous section), in this general case we have found only 
global minimum points of the TDP (|2"TT) of the form (Mo 7^ 0, Ao = 0) (the chiral symmetry breaking phase II) and 
(Mo = 0, Ao 7^ 0) (the superconducting phase III). Therefore, it is useful to present the analytical expressions for the 
thermodynamic quantities (|27|) in each of the phases II and III. Namely, for the phase II we have 



1 

I phase II 

_ MB 
I phase II 277- 

Ms 



- M$)0(in - M ) + sign( m )(^ 2 ~ M 2 )^(|/i 2 | - M )] 



- M o) 



I phase II ^ 

whereas for the phase III these quantities look like 
I 1 

I phase III 2 77 



- M ) - Bign{pa){A - M SM\^\ - M )] , 
\^\0(\^\ - M )} , 



(31) 
(32) 

(33) 



Mi\/« + A o + A o In 



Mi + VMi + A o 



+ • s '.'/"(/'2 - 1 H x/Ml + + A 2 In l^l + V^ + A^ 



(34) 



_ Ms 

I phase III 277 



M1A/M1 + Ag + A 2 In 



Mi + + A 2 
A 



si 5 n(/z 2 ) IM2K/M2 + A^ + A^ln 



Im 2 | + Vm1 + a^ 



A n 



mI 



phase III 



Mi + Aq + JmI + A 2 



(35) 
(36) 



The case g\ > 0. Our investigations show that in this case at arbitrary nonzero values of /x and B the supercon- 
ducting phase is realized in the system. Since at 52 < the phenomenon takes place even at \x = and B = 0, we can 
say that the chemical potential and in-plane magnetic field enhance superconductivity, which was originally generated 
in this case by a rather strong interaction in the fermion-fermion channel (Gi > G c ). In contrast, at gi > the system 
is in the symmetric phase if fi = and B = 0. However, arbitrary small nonzero values of [i and B induce in this 
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SUPERCONDUCTIVITY 
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FIG. 11. The case fi ^ 0: The (/x, B)-phase portrait of the 
model at arbitrary fixed g\ < as well as at gi = — 1.5|gi|. 



FIG. 12. 27ie case / 0: The (/i, B)-phase portrait of the 
model at arbitrary fixed g\ < as well as at gi — 0.5|gi|. 
Here Hi and II2 denote the chiral symmetry breaking 
phases with n = 0, m — and respectively. 
The notation III stands for the superconducting phase. 



case the superconductivity. In Fig. 9 and 10 the behavior of the superconducting gap Ao and such thermodynamic 
parameters of the model, as particle density n, magnetization m and magnetic susceptibility x vs B are presented at 
arbitrary fixed g\ > and gi = 0.5gi as well as at fixed chemical potential value, fj, = 0.5/g\. 

The case g± < 0. We have found that in this case the (//, B)-phase structure of the model is richer than in the 
case gi > 0. Indeed, supposing that gi = — 1.5|pi j , where g\ is arbitrary fixed and negative, it is easy to find the 
(/x, B)-phase portrait of the model drawn in Fig. 11. There in the chiral symmetry breaking phase both the particle 
density n and the magnetization m are equal to zero. Qualitatively, the similar B)-phase structure occurs for each 
coupling 32 from the interval 52 € (— k\gi\, — Iffil), where k ~ 3.08. However, if 32 > (recall, g\ < 0) or 52 < — k\di\, 
then the situation is changed qualitatively. In this case for the representative choice of coupling constants, 32 = . 5 1 <7i | 
(31 is arbitrary fixed and negative), the (ji, i?)-phase portrait of the model is presented in Fig. 12. It is clear from this 
figure that chiral symmetry breaking phase II is divided into two regions denoted as Hi and II2 . In the region Hi the 
quantities n and m are still equal to zero, whereas in the phase II2 both n^O and m ^ 0. On the boundary between 
chiral symmetry breaking phases Hi and II2 the relation \g\\fi + ^b\qi\B = 1 is valid. Moreover, we represent in Fig. 
13 and 14 the plots of gaps Mq and Ao as well as particle density n, magnetization m and magnetic susceptibility x 
as functions of external in-plane magnetic field B at g^ = 0.5|(?i| and /1 = 0.7/|<?i|. 

Finally, we would like to note that the TDP (|21[) is symmetric with respect to the transformation \i -f> v. So, such 
physical quantities as Mq, Ao, x remain intact, whereas n O m if the transformation fi ■<-> v is performed. We remark 
that all figures 9, 10, 13, and 14 show these physical quantities only as functions of v = [IbB at some fixed values 
of fi, i.e. its dependence on fi at fixed v is not presented in an explicite form in our paper. However, there is no 
need for a special numerical calculations in this direction, taking into account the above-mentioned symmetry under 
the permutation fi o v. Indeed, let us suppose that g\ > 0, gi = 0.5<?i and hbB = 0. 5/gi. In this case, in order to 
consider a behavior of the quantities Ao, n, m and x vs it is sufficient to imagine that the horizontal axis in Figs 9 
and 10 corresponds to a variable fi. Then the lines 1 and 2 of Fig. 9 will mean the plots of Ao and m vs \x, whereas 
the lines 1 and 2 of Fig. 10 will correspond to the plots of n and x vs respectively. In a similar way it is possible 
to extract from Figs 13 and 14 the information about behavior of Mo, Ao, n, m and x vs /i in the case of g\ < 0, 
32 = 0.5|3i| and [i B B = 0.7/|#i|. 

Hence, as it follows from the above consideration, at /i 7^ and sufficiently strong external in-plane magnetic 
field the superconducting (or Cooper pairing) phenomenon appears in the framework of the (2+l)-dimensional GN 
model (1) (of course, if G2 > 0). In particular, if G1.2 < G c , i.e. g± 2 > 0, then superconductivity is induced by 
infinitesimal values of fj, and/or By. Moreover, it is clear from our investigations that superconducting phase of the 
model is accompanied by a magnetization, or spin polarization, with positive valued susceptibility. It means that 
induced magnetic moment (which disappeares at By = 0) of the system and external magnetic field -By have the 
same direction, i.e. the superconducting state is a paramagnetic one. It is not a diamagnetic, as in conventional 
superconductivity. Note, the magnetic superconductivity phenomenon (which includes both the paramagnetic and 
ferromagnetic superconductivity) has a long history of investigations in condensed matter physics (see, e.g., in [31]). 
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FIG. 13. The case fi / 0: The gaps M and A vs B 
at arbitrary fixed gi < as well as at 172 = 0.5|gi| and 
fi = 0.7/[<7i|. Curves 1 and 2 are the plots of the dimen- 
sionless quantities |<?i|Mo and |<?i|Ao, respectively. Here 
Hb \gi\B c ^ 0.937. 



FIG. 14. The case fj, ^ 0: Particle density n, magneti- 
zation m and magnetic susceptibility x vs B at arbitrary 
fixed gi < as well as at gi = 0.5|gi| and /x = 0.7/|<?i|. 
Curves 1, 2 and 3 are the plots of the dimensionless 
quantities g^n, g\mj [is and \gi\xl A, respectively. Here 
iXB |ffi \B C « 0.937. 



VI. SUMMARY AND DISCUSSION 



In this paper we considered the (2+l)-dimensional GN-type model (1) with chiral and superconducting interaction 
channels under the influence of an (in-plane) external magnetic field. It is shown that in-plane magnetic field By catal- 
izes a creation of the paramagnetic superconductivity in the system. It means that if at Bn = the electromagnetic 
U(l) symmetry group was not broken, then at \B\\ | > B c , where a critical field B c can be even zero at some particular 

relations between coupling constants G\ t 2, there is a spontaneous breaking of U{1). In addition, if at Bn — the U{1) 
symmetry was broken spontaneously (due to a rather strong interaction in the fermion-fermion channel), then nonzero 
values of B\\ enhance superconductivity, i.e. the superconducting order parameter Ao is an increasing function vs 

|_Bj||. In addition, in-plane magnetic field induces also a nonzero spin polarization (paramagnetic magnetization) in 
the superconducting state. 

Various particular cases of the problem have been already considered earlier. Indeed, twenty years ago, it was found 
that at G2 = and G\ > an arbitrary weak perpendicular external magnetic field B± can not only enhance the 
chiral symmetry breaking (if Gi > G c ) but also induce spontaneous breaking of the chiral symmetry at Gi < G c (see, 
e.g., [13]). In contrast, recently it was established [2l[ that an application of an external parallel to the system 
plane (in-plane) magnetic field -By results in the restoration of chiral symmetry (recall, in this case G2 = 0). For 
an explanation of such a different reaction of the chiral symmetry on perpendicular and parallel magnetic fields one 
should remember that B± acts on the orbital angular momentum of electrons, whereas the in-plane field -By couples 
only to their spin. In the last case, due to the Zeeman effect, increase of the electron chemical potential takes place 
which eventually leads to the restoration of the chiral symmetry 0-0, EH . 

On the contrary, the results of our present paper demonstrate that the response of the electromagnetic U(l) sym- 
metry of (2+l)-dimensional GN-type models to the applied external magnetic field is completely different. Indeed, 
it was shown, in the framework of the (2+l)-dimensional GN-type model with nonzero coupling G2 of the super- 
conducting channel and at G\ = 0, that restoration of the U(l) symmetry takes place at sufficiently strong external 
perpendicular field B± [25j . Moreover, the main result of our paper is that if only an external in-plane magnetic field 
is included in the (2+l)-dimensional GN model (1) (at least with G2 > 0), then at sufficiently strong values of Bn 
the superconductivity phenomenon appears, i.e. the spontaneous breaking of the electromagnetic U(l) symmetry is 
originated. In particular, if an interaction in both channels of the GN model (1) is sufficiently weak, i.e. G1.2 < G c 
(or (71,2 > 0), then superconductivity is induced by infinitesimal values of Bn. 

To understand the last property, one should take into account that due to the Zeeman effect, i.e. due to the Zeeman 
;/-terms in Lagrangian (1), an external in-plane magnetic field -By induces (at /i = 0) Fermi surface for electrons in 
any planar system. However, in our model, since G2 > 0, there exists an attractive interaction between electrons just 
above the Fermi surface. So, the Fermi surface is unstable in favor of Bose-Einstein condensate of Cooper pairs, and 
a new (superconducting) ground state is formed in the system. 

Finally, we would like to recall about some QFT examples in which similar effects, although in a quite different 
physical contexts, are also observed. First of all, it is important to note that an ability of strong external magnetic 
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field to induce electromagnetic superconductivity of vacuum was recently established in the framework of QCD (see, 
e.g., the review [Hj]). Next, it is argued in some papers (see, e.g., [33|]) that ferromagnetism and superconductivity 
might coexist in dense quark matter, thus explaining super strong magnetic fields of magnetars. El 



Appendix A: Algebra of the 7-matrices in the case of SO(2,l) group 



The two-dimensional irreducible representation of the 3-dimensional Lorentz group SO(2,l) is realized by the fol- 
lowing 2x2 7-matrices: 

_ (I 0\ -1 . _ / i\ - 2 . ( II I 



7 =^3= { Q _! J, 1 =i°r={i J. 7 = ^ = ^ _i J < ( A1 ) 

acting on two-component Dirac spinors. 
They have the properties: 

Tr(rr) = w v ; [r,r]^~2ie^ a ; rr = -^^ a ia+9^, (A2) 

where g^ v = = diag(l, — 1, —1), 7 Q = g a pl^ , £ 012 = 1- There is also the relation: 

Tr(Yl v T) = -2ie^ a . (A3) 

Note that the definition of chiral symmetry is a bit unusual in three dimensions (spin is here a pseudoscalar rather 
than a (axial) vector). The reason is simply that there exists no other 2x2 matrix anticommuting with the Dirac 
matrices 7" which would allow the introduction of a 7 5 -matrix in the irreducible representation. The important 
concept of 'chiral' symmetries and their breakdown by mass terms can nevertheless be realized also in the framework 
of (2+l)-dimensional quantum field theories by considering a four-component reducible representation for Dirac fields. 
In this case the Dirac spinors ip have the following form: 

(A4) 

with ipi , ipz being two-component spinors. In the reducible four-dimensional spinor representation one deals with 
(4x4) 7-matrices: 7^ = diagtf 1 *, — 7 M ), where 7 M are given in (|A1[) . One can easily show, that (//, v — 0, 1, 2): 

rr(yv)= 4 5 ^ ; r y = a ^ +g ^. 

<r» v = \[-f,-f] = diag(-ie^ a %, -ie^ a %). (A5) 

In addition to the Dirac matrices 7^ (/i = 0, 1, 2) there exist two other matrices 7 3 , 7 s which anticommute with all 
7 M (/i = 0, 1, 2) and with themselves: 

7 3 = (?;Q,7 5 =7V7V-(?;- J )' ( A6 ) 

with / being the unit 2x2 matrix. Quite common is the following definition of the space parity transformation P in 
the four-component spinor formalism: 

P : iP{t,x,y)^ 1 5 1 1 i>(t,-x,y). (A7) 
Finally note that in terms of two-component spinors ipi , ip2 the parity transformation P (|A7[) looks like 

P : x, y) -> i^i) 2 {t, -x,y); ^{t, x, y) ^^{t, —x, y). (A8) 

Such a definition of the space parity transformation is commonly used in (2+l)-dimensional theories with four- 
component representation for Dirac spinors (see, e.g., in [2F 



4 In ferromagnetic superconductors the spin polarization, or magnetisation, is originated spontaneously, i.e. in this case there is a 
spontaneous breaking of the rotational invariance. In contrast, in our model (1) the magnetization of the superconducting state exists 
only in the presence of external magnetic field Bu , which destroys the rotational symmetry of the model explicitly. 
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Appendix B: The path integration over anticommutating fields 



Let us calculate the following path integral over anticommutating four-component Dirac spinor fields q(x), q{x): 



1 = J [dq][dq] exp (i J d 3 x[qDq - j(q T Cq) - j(qCq T ) 



(Bl) 



where we use the notations of Section ITT1 and D = j u id v — M + fry . Note in addition, the integral I is equal to the 
argument of the ln-function of the formula @ in the particular case with only one four-component Dirac spinor field 
(generalization to two Dirac spinors is trivial). Recall, there are general Gaussian path integrals [29j : 



| exp 1 1 



[dq] exp 



d 3 x 



d 3 j 



1 



q Aq + f] q 



(det(A)) 1/2 



exp 



\qAf + fjq T )= (det(A)) 1 / 2 cxp ( - J / d 3 x fjA^fj 



d 3 x 



1 — 1 -T 



(B2) 
(B3) 



where A is an antisymmetric operator in coordinate and spinor spaces, and r)(x), fj(x) are anticommutating spinor 
sources which also anticommutate with q and q. First, let us integrate in (|B1| over q-fields with the help of the 
relation (|B2I) supposing there that A = AC, qD — rj T , i.e. r\ — D T q T . Then 

1 = (det(AC)) 1/2 J[dq]exp(- % - J d 3 xq[AC + D{AC)- X D T ]f} . (B4) 

Second, the integration over g-fields in (|B4I) can be easily performed with the help of the formula (|B3[) . where one 
should put A = AC + D(AC)~ 1 D T and fj = 0. As a result, we have 



I = (det(AC)) 1/2 (det[AC + D(AC)- 1 D T ]) 1/2 = (det[A 2 C 2 + DC' 1 D T C]) 1/2 . 
Taking into account the relations {d v ) T — —d v and C^ 1 {^ V ) T C — —^f (v — 0, 1, 2), we obtain from (|B5 



I = (dct[- A 2 + £>+£>_]) 1/2 = (dot B) 
where D± — ^ v id u — M ± fry . Using the general relation det-B = exp(Trln_B), we get from (|B6 

d 3 p 



,1/2 



In/ = iTrln (B) = ^ J ^ ln(A^)) J d 3 x. 



(B5) 
(B6) 

(B7) 



(A more detailed consideration of operator traces is presented in Appendix A of the paper [30].) In this formula 
symbol Tr means the trace of an operator both in the coordinate and internal spaces. Moreover, \i(j>) (i = 1,2) in 
(|B7[) are two twice degenerated eigenvalues of the 4x4 Fourier transformation matrix B(p) of the operator B, i.e. 



XiMp) = M 2 -p 2 1 -p 2 2 - f i 2 +p 2 ~A 2 ± 2J-NPp\ - APpj + APpl + fp\ + fp\ 



(B8) 
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